Disappearance of equilibrium with thermal spin bath induced by non-Markov process 
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In most small system weakly interacting with a thermal bath will finally reach the thermal 

state with the temperature of the bath. We show that this intuitive picture is not always true by 
a spin star model where the non-Markov effect predominates in the whole dynamical process. The 
spin star system consists of a central spin homogeneously interacting with an ensemble of identical 
noninteracting spins. We find that the correlation time of the bath is infinite, which implies that 
the bath has a perfect memory, and that the dynamical evolution of the central spin must be non- 
Markovian. A direct consequence is that the final state of the central spin is not the thermal state 
equilibrium with the bath, but a steady state which depends on its initial state. 

PACS numbers: 03.65.Yz, 03.67.-a, 73.21. La 



I. INTRODUCTION 

A quantum system can never completely be isolated 
from its surrounding environment, and in many cases it 
must be treated as open system The investigations 
on quantum open systems have attracted more and more 
attentions and many relative theoretical approaches have 
been developed, such as the Markov approximation 
|1| , the quantum jump approach [4, 5] , the quantum state 
diffusion @, 0] , the canonical transformation [1, Q as well 
as the path integral method [l(| • 

It is a common sense that, when a small system in- 
teracts with a large thermal bath, it will finally reach 
the thermal state in equilibrium with the bath. How- 
ever, we will show that it is not always true by studying 
the reduced dynamics of the central spin which homo- 
geneously interacts with an ensemble of identical non- 
interacting spins. Such spin star configuration [llT[20l | 
can be realized in many solid-state quantum systems, for 
example, the semiconductor quantum d ot [ 2ll 426l| or the 
nitrogen- vacancy center in diamond (27l . |28|. 

Generally speaking, the environment interactin g w ith 
the small system is composed of multimode bosons [10| or 
spins }29l. l30j. For such an environment, the correlation 
time of the environment is much shorter than the charac- 
teristic time of the system. Especially, under the Markov 
approximation, the correlation function is assumed to be 
a S type of time interval [H-Q , which implies that the en- 
vironment's memory effect can be neglected during the 
time evolution. However, in our spin star configuration, 
the correlation function of the bath is found to be a peri- 
odic oscillation function without damping. The thermal 
bath then has a perfect memory and the reduced dy- 
namics of the central spin must be non-Markovian. As a 
result, the central spin can not reach the thermal state 
in equilibrium with the bath. 

In this paper, we study the dissipation of the central 
spin in a spin star svstem. Compared to the models pro- 
posed in Refs. [U, [IH, the free Hamiltonians of the cen- 
tral spin and the thermal bath, which do not commute 



with the interaction Hamiltonian, are taken into account. 
The dynamics of the central spin is then dramatically dif- 
ferent from that without the free terms. In our system, 
the whole Hilbert space can be decomposed into a direct 
sum of smaller subspaces due to the exchange symme- 
try of the bath spins, and then we are able to study the 
dynamics of the system exactly. We find that the final 
state of the central spin is a steady state which depends 
on its initial state, instead of the thermal equilibrium 
state as obtained under the Markov approximation [3lj |. 
The initial thermal equilibrium bath state is composed of 
a minority of the eigenstates, and most of the bath states 
remain unchanged during the time evolution, which leads 
to a perfect memory of the bath. This explanation is con- 
firmed by serval distinct discrete peaks in the spectrum 
diagram. 

The paper is organized as follows. In Sec. UH we set up 
our model as a spin star configuration and investigate the 
reduced dynamics of the central spin in detail. In Sec. lIIIl 
we find that the central spin can not reach the thermal 
state in equilibrium with the bath, and we further analyze 
the underlying physical mechanism with the help of fast 
Fourier transformation. In Sec. IIV1 we demonstrate that 
the traditional Markov approximation is not applicable 
by investigating the correlation time of the thermal bath 
as well as the mutual entropy between the central spin 
and the thermal bath. In Sec. El we draw the conclusions 
and make several remarks . 



II. THE MODEL AND THE REDUCED 
DYNAMICS OF THE CENTRAL SPIN 



A. The Hamiltonian 

We consider a spin star configuration as shown in 
Fig. [1] The system is composed of N + 1 localized spin 
1/2 particles. The central spin is immersed in a bath of 
N spins which are arranged in a circle. 
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In our consideration, the reduced density matrix of the 
central spin p s is a 2 x 2 matrix in the basis of states 
{It) j It)} w ith It) an d \l) being the spin up and spin 
down states respectively. The diagonal element (t \p | t) 
corresponds to the probability in its spin up state. In 
what follows, we will focus on the evolution of the proba- 
bility and investigate the underlying physical mechanism. 



FIG. 1. (Color online)A schematic diagram of spin star con- 
figuration. The central spin with the energy level spacing cjo 
homogeneously couples to the surrounding bath spins with 
the energy level spacing uj. The bath spins do not interact 
with each other. 



The Hamiltonian of the global system is 
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where <j z and a x are the Pauli operators, ujq and uj are 
the frequencies of the central spin and the bath spins re- 
spectively. The first two terms in the r.h.s. of Eq. ([T]) 
are the Hamiltonians for the central spin and the bath 
spins, the last term describes that the central spin ho- 
mogeneously couples to the bath spins with the coupling 
strength g. Throughout this paper, we set h= 1. 

Since the Hamiltonian (JT|) is invariant with respect to 
the exchange of arbitrary two spins in the bath, it is 
reasonable to define the collective operators: 



N (i) N (i) 

^ = E^ = E^- 



Then the Hamiltonian ([TJ becomes 



a = f 4°) + luj z + 
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The central spin thus couples to a collective bath angular 
momentum. The angular momentum of the bath is a 
changeable one and the total angular momentum j can 
take values as j = N/2, N/2 - 1, • • • , N/2 - [N/2], where 
the notation [a] denotes the maximum integer not larger 
than a. It is obvious that the total angular momentum of 
the thermal bath is a conserved quantity, i.e., [J, a] = 0. 

Utilizing the exchange symmetry of the bath, we inves- 
tigate the reduced dynamics of the central spin under the 
influence of the thermal bath. The state of the central 
spin is described by the reduced density matrix 



p b (t)=Tr B [p bB (t)}, 



(4) 



where Tr b is the partial trace over the degree of freedom 
of the thermal bath and p (t) is the density matrix of 
the global system including the central spin and the ther- 
mal bath at time t. The density matrix p SB {t) undergoes 
the unitary evolution 



p bB {t) = exv(-iat)p bB (0) exp{iat) 



(5) 



B. The initial state of the thermal bath 

We assume that the initial state of the global system 
can be factorized into an uncorrelated product state: 



p bB (0)=p b (0)®p»(0) 



(6) 



with p s (0) and p B (0) being the initial density matrices 
for the central spin and the thermal bath respectively. 

The initial state of the bath is a thermal equilibrium 
state 



exp(-/3ujj z ) 
Tr[exp(-/3 W J Z )] ' 



(7) 



where the z-component of the collective angular momen- 
tum J z is defined in Eq. and j3 = l/fc^Tis the inverse 
temperature of the bath, with ks being the Boltzman 
constant, which will be set to unit in the following. 

The Hilbert space of the bath is an N— fold tensor 
product of two dimensional spaces. A direct diagonal- 
ization of the Hamiltonian is practically impossible when 
the number of bath spins is large. Fortunately, the 
huge Hilbert space C B can be decomposed into a direct 



sum of smaller invariant subspaces C^,. . 



(8) 



Here, the 2j + 1 dimensional subspace Cj^,. is spanned by 
the states with the total angular momentum j(j < N/2). 
We introduce the index jj = 1 , 2 • • • to label the dif- 
ferent invariant subspaces with the same j, where the de- 
generacy is induced by the different possible orientations 
of the spins in the thermal bath. A simple calculation 
gives the degeneracy 
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The above analysis naturally motivates us to rewrite 
the initial density matrix of the thermal bath p B (0) in 
the direct summation of the invariant subspaces. To this 
end, we define a set of orthogonal bases Y)j, j, m) in the 
subspace Cf^ . . These states are the eigenstates of J 2 with 



eigenvalues j (j 



1) and of J z with eigenvalues m. Then, 
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FIG. 2. (Color online)The relative weight W(M) versus the 
excitation numbers M in different temperatures. The number 
of bath spins is N = 31. We use f3 = 1/fcsT to represent 
the temperature of the bath. The corresponding temperature 
gradually increases from (a) to (f). All parameters are in the 
unit of ui. 



p B (0) can be written as 



p B (0) = -exp(-^J 2 ) 



l 3 ~<i= a i 

= ~^Y Y Y ex P(-^ WTO )l73»i>w}(7j,i,i 

3 m=—j 7j = l 
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where the partition function Z is defined as 



3 ^= a j 

z = Tr [Y Y Y ex P(-^ m )l7i>i.w)(7j,i,m|]. 

3 m=-j 7j=l 

(11) 

For the thermal equilibrium state given by Eq. (ITOl) . 
the relative weight for M excitations among N bath spins 
(i.e., M spins in their up states and the other N—M spins 
in their down states) in the thermal bath is 

WW-(») -*-»*f- 1 '™ (12, 

In Fig. [2j we plot the relative weight W(M)(N = 31) in 
different temperatures. It shows that, the states with low 
excitations play a predominate role in the low tempera- 
tures. As the rising of the temperature, the states with 
high excitations gradually emerge in the initial state, and 
the weights of states with low excitations are suppressed. 
When the temperature is high enough, the distribution 
is temperature independent (shown in Fig.f2](e) and (f)). 



We point out that, the excitation number is not a con- 
served quantity in our system. Therefore, the states with 
different excitations may couple to each other, generally 
depends on whether they belong to the same invariant 
subspace Cf y . . 



C. The reduced dynamics of the central spin 

In this subsection, we will study the reduced dynamics 
of the central spin under the influence of the thermal 
bath. We illustrate the dissipation of the central spin by 
calculating the probability in its spin up state. 

As discussed in the previous subsection, the Hilbert 
space of the thermal bath is composed by a direct sum 
of smaller invariant subspaces. Furthermore, the Hilbert 
space C SB of the global system (including the central 
spin and the thermal bath) is given by C SB = (Bj 1:j Cj B , 
where 



(13) 



with Cf 7j being defined in Eq. (JSJ and C s being the 2- 
dimensional Hilbert space for the central spin. In the 
basis of 7j, j, m) = ® 7^, j,m), the Hamiltonian is 

JV 

written as H = J2j J2^=i ^-n witn 
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where the matrix element ^™ is defined as 

= (t*';73,j,rn'\H\n;73,j,rn). (15) 

Correspondingly, we rewrite the initial den- 
sity matrix of the global system as p SB (0) — 

N 

EjE^i^p^yz with 
1 3 

(0) = -j— p s (0)(g) exp(-/3wm)|7 i ,j,m)(7 7 -,j',m| 

7n——j 

(16) 



P SB 



31j 



where 

3 

z n =Tr [ Y ex P(-Pum)\'y j J,m)('Yj > j,rn\] (17) 

m=—j 

is introduced to guarantee the unity trace of Pj B (0). In 
the above equations, 'H.j lj and Pj B (0) can be understood 
as the projection of the Hamiltonian and the initial den- 
sity matrix on the subspace respectively. 

In the subspace Cj B , the "density matrix" pj B expe- 
riences the unitary evolution 

ijljtjPjr/j 



p? B (t) = exp(-iH n t)p s , B (0) exp(^ 7 . t) . (18) 
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and the reduced density matrix for the central spin is 
Pj 7:i (t) = T^Bj^iPjj^t)], where Ti Bj7j [•] represents the 
partial trace over the degree of freedom of the thermal 
bath in the subspace C^. 

Hence the probability for the central spin in its spin 
up state is 



Pit) 



N ry 
a j Z 



^(tl^WIt). 



(19) 



III. RESULTS AND DISCUSSIONS 

In general cases, a small system interacting with a large 
thermal bath will be thermalized. In other words, the 
small system will reach the thermal state in equilibrium 
with the bath. For example, for a single spin-1/2 particle 
as the small system which is described by the Hamilto- 
nian H = wo<r z /2, the final density matrix would be 



p(t 



exp(-^w er z /2) 



' Tr[exp(-/3w <7 2 /2)] 
and then the probability in its up state is 
p (t = oo)= l_tanh(^ /2)_ 



(20) 



(21) 



It shows that the probability depends only on the temper- 
ature of the thermal bath but not on the initial state of 
the spin. Such results can be obtained under the Markov 
approximation [3lLl3^. In the following, however, we will 
show that this intuitive picture is not true in our system. 

Through a straightforward numerical calculation based 
on Eq. (fT9|) . we illustrate the dynamics of the probability 
for the central spin in its spin up state in Fig. [3] In the 
figure, the red (blue) curve represents the case when the 
central spin is in its up (down) state initially. It shows 
that, the final state of the central spin depends on its 
initial state. As shown in the figure, the probability os- 
cillates around a fixed value and the amplitudes of the 
oscillation decrease with the increase of the bath tem- 
perature. In other words, the probability tends to be a 
constant when the temperature of the spin bath is high 
enough. 

Furthermore, when the central spin initially prepared 
in the spin up or spin down state, we will have (<Tx) = 
(<jy ) = during the evolution. This can be clarified 
in the viewpoint of Zi symmetry of the system. By the 
analogy of Ref. |33j . we define the parity operator V = 
exp[iTr(J z + ai°^ /2)] which is conserved with respect to 
the Hamiltonian (J3J, namely, [V, H] = 0. In our consid- 
eration, the initial state has a fixed parity and it will stay 
in the same parity during time evolution due to the con- 
servation of the parity. However, the central spin opera- 
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FIG. 3. (Color online) The probability for the central spin in 
its spin up state as a function of time t. The initial state of the 
central spin is spin up state for red lines and spin down state 
for blue lines. The parameters are set as the loq = l,g = 0.1. 
The inverse of temperature is (a)/3 = oo, (b)/3 = 3, (c)/3 = 1, 
(d)/3 = 0.5, (e)/3 = 0.1, (f)/3 = 0. The number of spins in the 
thermal bath is N — 31. All parameters are in the unit of u). 



Therefore, their average values will keep zero at any time 
t. 

The above discussions show that, when the thermal 
bath is large enough, the central spin will evolve into 
a state which satisfies (<Jx) — = and (eri°') ~ 

constant. Therefore, we can safely conclude that the cen- 
tral spin will reach a steady state instead of the thermal 
equilibrium state. 

As demonstrated above, the central spin can not be 
perfectly thermalized as expected by interacting with an 
ensemble of identical spins. The main reason is that, 
the initial thermal state of the bath is composed of a 
minority eigen states of the Hamiltonian, and therefore 
the energy level transitions during time evolution only 
occur among the few eigen states which are related to 
the initial state, leaving other states unchanged. As a 
result, the bath has a perfect memory and the initial 
thermal equilibrium is broken as the time elapses, which 
leads to the disappearance of the equilibrium state of the 
central spin. 

To investigate the energy level transition in detail, we 
switch, by the fast Fourier transformation (FFT) tech- 
nology, from the time domain to the frequency domain, 
where we can obtain the transition frequency quantita- 
tively. 

We plot the frequency spectrum diagram for the prob- 
ability obtained in Eq. (|19l) in different temperatures in 
Fig. [4] assuming that the central spin is in its spin up 
state initially. It follows from Fig. @] that, there exist 
two distinct peaks in zero temperature(/? = oo). As the 
temperature rises, the high frequency peak gradually de- 
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FIG. 4. (Color online) The frequency spectrum of the proba- 
bility in Eq. (|19[) by FFT. The parameters are same as that 
in Fig. |31 The central spin is in the spin up state initially. 



creases until disappears. Meanwhile, more and more low 
frequency peaks with relative weak strength emerge. Ac- 
tually, all of these peaks correspond to the related energy 
level transitions. 

In the low temperatures, nearly all of the bath spins 
are in their spin down states initially. The initial state 
only locates in the subspaces Cf 7j with large j, for ex- 
ample, j = N/2, N/2 - 1, N/2 - 2. Moreover, in the case 
of weak coupling, the transitions only occur among a few 
lowest intra-subspace energy levels. This is similar to 
the Rabi model, which describes the interaction between 
a two- level system and a single mode bosonic field [34j . 
In zero temperature, the Rabi oscillation like phenomena 
recurs in Fig. EJa) where we note that the "anti-rotating 
wave" term leads to the inequality of the oscillation am- 
plitudes. As for the case of the low temperatures, it also 
reproduces the phenomena of collapse and revival, which 
can be observed clearly in Fig. Elb). 

As the temperature further rises, the components of 
initial state emerge in more and more subspaces C^, 
with small j and the relative weights in the subspaces 
with large j gradually decrease (shown in Fig.0}. In these 
small j subspaces, the energy level differences are smaller 
than those in the large j subspaces. Therefore, the peaks 
with high frequency gradually disappear and the low fre- 
quency peaks arise(shown in Fig. When the temper- 
ature is high enough(for example when j3 < 0.1), the 
frequency spectrum and the dynamical evolution process 
are temperature independent. We see that the curves for 
= 0.1 and = coincide with each other in Fig. 2] and 
Figs. He, f). 

The frequency spectrum in Fig. 2] shows that the prob- 
ability (the r.h.s. of Eq. (TT9l) ) is a superposition of oscil- 
lation functions with different frequencies and all of the 
components evolute with time in unison during a short 
time. However, the interfere effect among different com- 
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FIG. 5. (Color online)The fluctuation of the probability AP 
as a function of the bath spins' number N for different tem- 
perature. The parameters and initial state are same as that 
in Fig. 4. The case of = 0.1 and = completely coincide 
with each other. 

ponents becomes important once the time is long enough, 
and the probability oscillates around a fixed value as the 
time elapse. The different behaviors of the probability in 
two time scales are clearly shown in high temperatures 
in Fig.H (d-f). 

To depict the quantum interference phenomena quan- 
titatively, we define the probability fluctuation 



AP = (p(i) - P(t)) , (22) 

where a denotes the mean value of a over time. In Fig.[5j 
we plot the curve of AP as a function of number of bath 
spins in different temperatures. It shows that AP de- 
creases monotonously as the increase of the spins number 
and the temperature of the bath. Both of the increasings 
imply that more components participate the summation 
in Eq. (fl"9|) . it is the destructive interference of the dif- 
ferent components that weakens down the fluctuations. 

IV. THE INVALIDITY OF THE MARKOV 
APPROXIMATION 

In the study of the open system, the Markov approx- 
imation is widely used to discuss the reduced dynam- 
ics of the system. The Markov approximation usually 
works well when the small system weakly couples to a 
large thermal bath. Under the Markov approximation, 
the typical time of the internal correlation of the bath is 
assumed to be much shorter than the characteristic re- 
laxation time of the system, which implies that the bath 
correlation decays more rapidly compared to the system 
damping. As a result, the bath loses its memory and does 
not change significantly with time. Therefore, the system 
does not affect the bath significantly, but the bath does 
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affect the system. However, the Markov approximation 
loses its effectiveness in our system even if the conditions 
of large thermal bath and weak coupling are both sat- 
isfied. In this section we will show the invalidity of the 
Markov approximation in two aspects: the internal corre- 
lation time of the bath tends to be infinite and the central 
spin does correlate with the bath all the time during the 
time evolution. 

Firstly, we study the correlation of the thermal bath. 
It is convenient to work in the interaction picture with 
H = Wo<?z /2 + w 5^02/2, the internal correlation 
function of the thermal bath is calculated as 

2 cosh(pcj/2) 

with T(t) — gJ2j <T -' >e ~ l ^ jt an d (-)b being the average 
over the thermal equilibrium state of the bath in temper- 
ature T which is described by the density matrix 



exp 



(-/fc>5>w/: 



Tr exp(-/3w/2^a^/2 



(24) 



Here, Tr means the trace over the degree of freedom of 
the spin bath. 

It is shown that the correlation function in Eq. (|23|) 
is a periodic oscillation function without any damping. 
Therefore, the memory effect of the thermal bath can 
not be neglected and the condition of the Markov ap- 
proximation is naturally broken. The same conclusion in 
spin star system is also obtained by a direct analysis of 
the Redfield master equation in Ref. fTTj j . 

Secondly, we will discuss the correlation between the 
central spin and the thermal bath applying the concept of 
the mutual entropy. In quantum information, the mutual 
entropy between two subsystems (the central spin and the 
thermal bath in our system) is defined as: 



s( P b ) + s{ P ») - sgo 



(25) 



where S = — J2i ^« l°g2 C^i being the eigenvalues of p) 
is the von Neumann entropy (35l . l3o | . Under the Markov 
approximation, the influence of the open system to the 
bath is ignored, and the global system will stay at the 
tensor product state 



p bB {t)=p b {t)®p? h 



(26) 



with pf h being the initial thermal equilibrium state of 
the bath. The mutual entropy of such product state is 
zero because the correlation between two subsystems dis- 
appears. However, this is not the case for our system. 
In Fig. [5J we plot the mutual entropy obtained by our 
approach in different temperatures. It shows that the 
mutual entropy never disappears, which implies that the 
central spin does correlate with the thermal bath, and is 
incompatible with the result from the Markov approxi- 
mation. This also verifies the invalidity of the Markov 
approximation. 
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FIG. 6. (Color online)The mutual entropy between the cen- 
tral spin and the thermal bath as a function of time at dif- 
ferent temperatures. The parameters and the initial state are 
same as that in Fig. [3] 



V. CONCLUSIONS AND REMARKS 

In this paper, we have studied the reduced dynamics of 
the central spin in a spin star system and calculated the 
probability of the central spin in its spin up state. We 
find that the probability oscillates around a fixed value 
and the oscillation amplitudes decrease with the increase 
of temperature and the number of bath spins. To inves- 
tigate the energy level transitions during the time evolu- 
tion, we utilize the FFT technique to obtain the spectrum 
diagram. In the low temperatures, there only exist two 
distinct peaks and the probability function shows peri- 
odicity which is similar to the Rabi oscillation. As the 
increase of the temperature, it exhibits more peaks in 
the spectrum diagram, which shows that the probability 
function is composed of finite components with differ- 
ent frequencies. It is the destructive interfere of different 
compositions that weakens down the oscillation. 

Here we emphasize that the reduced state of the cen- 
tral spin can be constructed from its probability in the 
up state. This is because the average values of the x- 
component and y-component keep zero during the time 
evolution, which is directly proved from a consideration 
of parity conversation. 

In our system, the conventional Markov approximation 
is not reasonable even in the case of weak system-bath 
coupling. On one hand, it is discovered that the correla- 
tion time of the bath tends to be infinite, which implies 
that the bath has a perfect memory. On the other hand, 
the correlation between the central spin and the thermal 
bath also implies the invalidity of the Markov approxi- 
mation in our system. Due to the non-Markov effect, the 
central spin can not be perfectly thcrmalizcd by interact- 
ing with a bath of identical spins. Instead, it will evolute 
into a steady state which depends on its initial state. 

Technically, we point out that the Hilbert space is too 
huge to make a direct numerical calculation in our system 
when the number of bath spins is very large. However, 
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the huge Hilbert space can be decomposed into a direct 
sum of smaller subspaces utilizing the exchange symme- 
try of the bath spins and we are able to study the dynam- 
ics of the system exactly. By using of the same approach, 
the entanglement and other problems in quantum infor- 
mation related to spin star system can be solved without 
any approximation or specific numerical technique. 

In summary we give a spin star system to demonstrate 
a quantum open system without the Markov approxima- 
tion for whatever weak coupling between the system and 
the environment. In our system, a considerable correla- 
tion between the system and the environment exists even 
in a long time, and the central spin will finally evolve 
into a steady state dependent of its initial state. We 



hope that our investigation on this simple system will in- 
crease our understandings on the dynamics of quantum 
open system, a central topic in the physical realizations 
of quantum information processes. 
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